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A general recipe to define, via Nother theorem, the Hamiltonian in any 
natural field theory is suggested. It is based on a Regge-Teitelboim-like 
approach applied to the variation of Nother conserved quantities. The 
^ik ' Hamiltonian for General Relativity in presence of non-orthogonal bound- 

h»«, ' aries is analysed and the energy is defined as the on-shell value of the 

("^ , Hamiltonian. The role played by boundary conditions in the formalism is 

t*~« ' outlined and the quasilocal internal energy is defined by imposing metric 

c 2 ^ I Dirichlet boundary conditions. A (conditioned) agreement with previ- 

ous definitions is proved. A correspondence with Brown- York original 
formulation of the first principle of black hole thermodynamics is finally 
^ , established. 

tu;)! 1 Introduction 

One of the most popular approaches to black hole thermodynamics is due to 
S^ ■ Brown and York; see ||l|, ||. They defined, through an Hamilton-Jacobi analysis 

JH I of a gravitational system, the quasilocal energy, the quasilocal angular momen- 

tum and the spatial stress of a spatially bounded region of spacetime. Those 
quasilocal surface densities were subsequently involved in describing the first 
law of black hole thermodynamics through a statistical analysis based on path 
integral methods. As a result, in the final expression of the first law, quasilocal 
physical quantities are involved rather than the global ones calculated asymp- 
totically. Some advantages may, in fact, be gained by considering spatially 
bounded self-gravitating systems instead of considering the properties of the 
whole universe. From a mathematical viewpoint, if we analyze only a spatially 
bounded region, the asymptotic behavior of the gravitational field becomes irrel- 
evant and different gravitational solutions with different asymptotic behaviors 
can be handled on equal footing. On the other hand, from a physical viewpoint 
it was noticed that a system in thermal equilibrium must feature a finite spatial 
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extent: a system of infinite spatial extent at fixed temperature is thermodynam- 
ically unstable (see e.g. [^). Hence, one could argue that, in establishing the 
first principle of black hole thermodynamics, the role of internal energy should 
be played by a suitable defined "surface energy" rather than by the total mass 
(calculated in asymptotic regions). Correspondingly, if the temperature is iden- 
tified with the intensive variable which is thermodynamically conjugated to the 
energy, it has to depend somehow on the surface rather than to be identified 
with a constant parameter. The same reasoning can be applied to all the rele- 
vant physical quantities which are expected to appear in the first law of black 
hole thermodynamics. This is the idea at the basis of the aforementioned papers 

0,1. 

Nevertheless the formalism developed in those papers deals only with folia- 
tions of spacetime with spatial boundaries orthogonal to the timelike boundary. 

In |4| |5| the assumption of orthogonal boundaries was removed and the 
quasilocal formalism of Brown- York was extended in order to deal with the 
general case. The formalism there developed starts from the so-called Trace- 
K gravitational action functional as introduced by Hayward in |^. Roughly 
speaking, the Trace-K action differs from the standard Hilbert action by metric 
dependent boundary and corner terms and it features a variational principle 
with fixed boundary metric. The canonical analysis of H, M goes on through an 
ADM splitting of the Trace-K action functional. The Hamiltonian ensues from 
the canonical reduction of the action; the energy (or it would be better to say 
the mass) is defined as the on-shell value of the Hamiltonian and turns out to be 
a pure boundary term. By focusing just on the (space -I- time) foliation of the 
timelike boundary, the formalism turns out to be independent on the foliation 
of spacetime inside or outside the boundary. The expression of the quasilocal 
energy for foliations with non-orthogonal boundaries was applied in Q to study 
the Lorentz-like transformation rules of the quasilocal energy with respect to 
boosted observers. In particular, the energy measured by observers infalling 
towards a black hole was calculated. Moreover, the close analogy between the 
quasilocal formalism and the thin shell formalism of Israel (see 0) was analysed 
in [y with the aim of physically supporting and motivating the interpretation 
of energy attributed to the mathematical formula for energy there suggested. 

A similar formulation was also given by Hawking and Hunter in [g| . However 
the analysis there developed was carried on by considering the foliation of the 
spacetime as a whole rather than focusing just on the foliation of the bound- 
ary. The Hawking-Hunter definition of the gravitational Hamiltonian features 
a "tilting term" which has the disadvantage to depend on the intersection angle 
of the boundaries. In order to remove this unaccettable term a careful choice of 
the background has to be imposed; see p|. 

Another definition of energy for non-orthogonal boundaries can be also found 
in Kijowski's paper M . Such a definition was there obtained by means of a sym- 
plectic analysis of the gravitational Lagrangian and through boundary Legendre 
transformations. The main contribution of Kijowski's approach is related to the 
way in which different boundary conditions are related to different definitions of 
energy (i.e. internal energy or free energy). Starting from the same gravitational 



Lagrangian, namely the Hilbert Lagrangian, different definitions of energy are 
recovered in through different control modes of the boundary data. Moreover, 
different boundary conditions reverberate in different Legendre transformations 
or, in other words, in different boundary symplectic structures. 

Many attempts to describe the quasilocal stress-energy content of a spatially 
bounded region of spacetime via Nother theorem can be also found elsewhere; 
see, e.g. [^ |Il| |l^, |T^ and references quoted therein. Nevertheless it was long 
but wrongly believed - see Q ~ that a Nother framework for conserved quanti- 
ties is not well-suited to deal with foliations of spacetime with non-orthogonal 
boundaries. The purpose of this paper is to contradict this common belief. In- 
deed, we shall show that a Regge-Teitelboim-like approach to Nother theorem 
allows to handle boundary terms in such a way that quasilocal Nother charges 
can be defined. Roughly speaking, the leading idea is to define the variation 
of conserved quantities 6Q along a 1-parameter family of field configurations 
rather than defining the conserved quantities Q in their own. It is in fact well- 
known that absolute conserved quantities do not have a proper physical meaning 
while it is instead meaningful to define relative conserved quantities, i.e. con- 
served quantities with respect to a fixed solution chosen as a background and 
which is introduced with the aim to define the "zero level" for conserved quan- 
tities themselves; see [^ 14, |l^, |l^, 0. The background solution has to be 



chosen carefully (and differently case by case) because it has to satisfy appro- 
priate boundary matching conditions. However it is also well-known that the 
problem of finding a suitable background, if possible, can be rather cumbersome 
in direct applications. Nevertheless, at a theoretical level, the problem can be 
bypassed by defining the variation SQ as suggested above. We remind that to 
each infinitesimal Lagrangian symmetry ^ one can associate via Nother theorem 
a conserved quantity Q{^) which is made of a bulk term Q|3^i]^(C) plus a surface 
term QsmiaceiO- O'^ ^ spacctime M of arbitrary dimension m (m ^ 2) the bulk 
term is obtained by integrating on a {m — 1) dimensional region E a suitable 
(m — l)-form £{^) which is called the reduced Nother current and which is van- 
ishing on shell, i.e. by setting Q^^i^(C) = I^^iO- The surface term Qsurfaco(0 
is instead obtained by integrating on the boundary dT, the (to — 2)-form U{£,) 
which is called the superpotential and satisfies a canonical splitting £ = £ + dU, 
where £ is the generator of Nother currents (see e.g. |l^, |l9[), i.e. by setting 
Qsurfacc(0 = Igs^iO- The supcrpotential is algorithmically calculated for ev- 
ery natural or gauge-natural theory starting from the Lagrangian of the theory 
and the vector field ^; see jlj, Q. If we consider a Cauchy hypersurface E 
in spacetime and a vector field ^ transverse to E the Hamiltonian can be a 
priori defined for a Lagrangian theory as the Nother charge relative to ^, i.e. 
H{(,) — Q{^)- Namely, the Hamiltonian is the Nother charge associated to the 
"time" desplacement vector field ^, having identified the (local) flow parameter 
of ^ with the (local) time. Nevertheless this definition does not produce the 
expected values when specialized to specific solutions (it is indeed affected by 
anomalous factors or divergence problems) mainly because of the fact that such 
a definition does not take into proper account a reference background. There- 
fore, this prescription leads to a definition of conserved quantities which has to 



be somehow corrected. To this end one defines the variation 5Q of the corrected 
conserved quantity Q by adding to SQ a suitable additional boundary term. 
The latter ensues from the variation 5Q of the bulk term and it keeps into 
account the behavior of the fields at the boundary. In this way the formalism 
resembles the original idea of Regge-Teilteiboim (see p^|), i.e. boundary terms 
arising in the variation of the bulk Hamiltonian are added (with a minus sign) to 
the very definition of Hamiltonian so that the variation of the new Hamiltonian 
does not contain boundary terms at all. However, once the variation 6Q of the 
corrected conserved quantity is defined, it remains to analyse the problem if 
the variation 6Q can be formally integrated. We shall see that this problem is 
tightly related with the boundary conditions we choose (i.e. it is related to the 
control boundary data). When SQ is integrable we obtain the conserved quan- 
tity Q — Qo defined up to a constant of integration Qq. The latter can be fixed 
as a zero level for the conserved quantity or, in other words, as a background 
reference. Namely, Qo corresponds to the Nother charge relative to a solution 
inside the set of solutions satisfying the same boundary (or asymptotic) condi- 
tions (and, in turn, it fixes the boundary or asymptotic conditions themselves). 
Notice, however, that the suitable matching conditions between a solution and 
its background are required only on the boundary (or at infinity) so that the 
background is allowed to have a different topology (for example, this is the case 
if we fix the Minkowski background for the Schwarzschild solution; see |13| for 
details). 

The variation of Nother conserved quantities was already considered in |^ . 
This approach was at the basis of the so-called covariant ADM formalism. The 
formalism has been deeply analysed and tested in its applications to General 
Relativity and it was shown that it reproduces the expected physical values of 
conserved quantities. Nevertheless, as far as we know, all the examples so far 
analysed in General Relativity cover only the non-orthogonal boundary case. 
It is our goal in this paper to generalise and extend the aforementioned for- 
malism in order to fill this gap, also in view of a unifying perspective on the 
various formalisms recalled above and sometimes apparently unrelated. To this 
end a new additional surface term has to be added to the Nother-based def- 
inition of Hamiltonian when dealing with non orthogonal boundaries. Indeed 
the variation of the corrected Hamiltonian H{(,) = Q(C) gives rise to an ad- 
ditional boundary term which is instead vanishing in presence of orthogonal 
intersections. Again, this term has to be added to H with a minus sign so 
that, eventually, the variation of the Hamiltonian does not really contain sur- 
face terms. It is also remarkable that the final formula we shall obtain for the 
Hamiltonian is independent on divergence terms added to the Lagrangian, so 
that it depends only on the field equations content of the equivalence class [L] 
of all such dynamically equivalent Lagrangians. From a mathematical view- 
point this is clearly an advantage. Indeed, we do not have to care about surface 
terms in the action functional. On the contrary, other Hamiltonian definitions 
based on a framework different from the Nother one (such as the ones based on 
the canonical reduction of the action functional or based on a Hamilton- Jacobi 
analysis) are deeply sensitive to boundary terms. In those formulations, in fact. 



boundary terms reflect into boundary conditions which have to be satisfied by 
the dynamical fields. Different surface terms in the action functional then cor- 
respond to different boundary conditions and, in the end, to different definitions 
of energy. Indeed, it is important to notice that for most physical systems there 
may exist different kinds of energy, each one corresponding to different choiches 
of boundary conditions and/or of control variables. For istance, if we agree 
that the internal energy of a system is the energy corresponding to Dirichlet 
boundary conditions, when we want to calculate the internal energy for a solu- 
tion of Lagrangian field equations, we have to go back to the Lagrangian itself 
and to look for divergence terms which have to be added to the Lagrangian to 
obtain the appropriate Dirichlet boundary conditions (in General Relativity the 
sought-for Lagrangian corresponds to the Hayward Lagrangian ||]). Moreover, 
when we try to consider a self-gravitating system as a thermodynamical sys- 
tem, different boundary conditions reverberate on different choices between the 
fields which have to be considered as intensive or extensive thermodynamical 
variables; see [0. Accordingly to the choice made on thermodynamical variables 
(i.e. on the boundary control modes of the system) one has then to exibit differ- 
ent action functionals, e.g. the microcanonical action functional, the canonical 
partition function or the grand canonical partition function. In this way the 
initial Lagrangian must be modified each time according to the model one is 
dealing with. 

On the other hand, the definition of the Hamiltonian based on Nother the- 
orem is unique inside each class [L] of Lagrangians, where the elements of each 
class [L] differ from each other only by total divergences^. No ad hoc physical or 
mathematical prescription is then required to select a well-defined representa- 
tive inside [L] because all of them lead to the same Hamiltonian. Nevertheless, 
if the Hamiltonian is unique for any element of the class, the question can arise 
how different definitions of energy can be recovered in this framework. In other 
words, how different boundary conditions can be imposed on the same Hamilto- 
nian? (A similar problem was analysed by Kijowski in [^ and also in ||l^, [l5||). 
In our case the answer has to be looked for into the way we consider the variation 
of the Hamiltonian. Namely, we can consider different one-parameter families of 
solutions each one containing the same solution (i.e. in the infinite dimensional 
space of solutions we consider different curves passing through the same point). 
Different one-parameter families of solutions (i.e. different curves) correspond 
to different ways in which the variation of the Hamiltonian is carried out and, 
in turn, to different definitions of energy. 

The plan of the paper is as follows. In Section Q we shall present the theoreti- 
cal framework leading to the definition of Hamiltonian and energy in relativistic 
field theories. We stress that even though the formalism will be presented with 
the aim of being applicated to General Relativity, it will be developed inside a 
geometrical framework which is common to all natural as well as to all gauge- 
natural relativistic theories so that, at least in principle, all classical theories of 
physical interest can be covered. Starting from Section we shall focus only on 



^ This is a cohomology class in the proper variational cohomology theory; see, e.g. |22l BS 



General Relativity which for simplicity will be considered only in vacuum. The 
generalization to cases with matter is an easy task which adds nothing to the 
understanding of the problem and, as such, we leave it to some interested reader. 
It is in fact the second order Hilbert Lagrangian (as well as any Lagrangian of 
order fc > 1) which is sensitive to the formalism presented here; first order mat- 
ter Lagrangians can in fact be treated in the usual way; see |14, ^. First of 



all we shall shortly review the geometry of the non-orthogonal foliati on o f (a 



portion of) spacetime and we shall fix the notation. Then, in subsection 3.1, the 
variation of the Hamiltonian together with Hamilton equations are calculated 
in order to test the formalism. The quasilocal energy of the gravitational field. 



defined as the on-shell value of the Hamiltonian, is analysed in subsection 3.2 
It will be there shown that the variation of the energy can be formally integrated 
provided Dirichlet boundary conditions are imposed. In this way a definition 
of quasilocal internal energy for a spatially bounded gravitating system is ob- 
tained whence the reference background subtraction terms are properly taken 
into account. Moreover, a comparison with the alternative definitions of inter- 
nal energy given in |j, g, |^ is carri ed o ut and a (conditioned) correspondence 



is provided. Finally, in subsection 3.3, the first principle of black-hole ther 



modynamics is easily obtained as a trivial consequence of the cohomological 
properties of Nother charges. 

2 Nother Theorem 

We shall here consider a field theory described by a Lagrangian of order k. 
According to the geometric approach to field theory (see for instance ||l^, ^, 
p5| , ^ ) we are assuming that the configuration bundle of the theory is a 
bundle (Y,M, tt). The base manifold M, with dimAf = m, is the space of 
parameters of the theory. From now we shall identify the manifold M with 
the physical spacetime. Fibered adapted coordinates on Y are {x'^,y^), where 
H = 1, . . . ,m = dim M, i — 1, . . . , n = dim Y — dim M . A configuration is a 
(local) section a : M — > Y, locally x^ i-^ {x^ ,y^ — a'^{x^)), which to each point 
X G M associates the values of the fields "y" (of components y*). Thence, a 
configuration completely describes the evolution of the fields which represent a 
given physical model. The fibered manifold {J'^Y, ir^ ,M) denotes the k jet pro- 
longation of the configuration bundle (i.e. the bundle where fields live together 
with their spacetime derivatives up to order k included) and j'^a : M — > J'^Y 
denotes the jet prolongation of the section a (i.e. j'^ denotes derivation up to 
order k). Natural fibered coordinates on J''Y are obviously defined and de- 
noted by {x^,y^, . . . , ?/^^...^j^)- In the framework so far developed a Lagrangian 
of order k can be seen as an horizontal form (or equivalently as a base pre- 
serving fibered morphism) L : J^Y — i-A™(M). It can be locally written as 
L = C{y\ . . . , yli-^,,,fj_^) ds where £ is a scalar density, called the Lagrangian den- 
sity, and ds = dx^ A . . . /\ dx"^ is the local volume element of m-dimensional 
spacetime. Once the Lagrangian is evaluated on (the jet prolongation of ) a 
section it gives rise to a m-form {j^a)*L on spacetime M . 



A vertical vector field X is a section of the vertical bundle VY — > Y, i.e. it 
is a vector field on the configuration bundle Y which is everywhere tangent to 
the fibers (or, equivalently, projects down to zero under the tangent map of the 
bundle projection). Locally, it reads as X = X^d/dy^ — [5y^)d/dy^ and it de- 
scribes what is classically called the "variation" Sy"^ of the dynamical fields. The 
fc-jet prolongation J^X is a section of the bundle J^{VY) — >Y and, roughly 
speaking, it describes the variation dy^,6y^^, . . . 5y^^_^ ^^ of the dynamical fields 
together with their derivatives up to order k. 

Given a vertical vector field X and a section a and denoting by (pT the 
(local) fiow of X we define a one-parameter family Gr of sections through the 
rule Gr = 4>T ° CT- The variation 6x^ with respect to X of a differential form lj 
on any J^Y [h > 1) is defined according to the rule 

Sxuj = -7-uj{j'^ar) 

From now on, for notational simplicity, we shall assume that the Lagrangian, 
as well as all horizontal forms on any prolongation of the configuration bundle, 
are already computed along sections so they actually define differential forms 
on spacetime, e.g. L stands shortly for {j^(j)*L. In other words we use the 
shorcut y* to actually denote y* = cr*(a;), i.e. the i-th component of the section 
a with respect to the i-th fibered coordinate y* . Analogously y* is used instead 
of di_i(T^{x) and so on for the derivatives j^y of higher order. Nevertheless we 
emphasize that all calculations we shall carry on in the standard language of 
differential calculus on the base manifold M can be alternatively carried on at 
the level of jet bundle calculus and only at the end the results so obtained can be 
pulled back on spacetime. In the latter case globality and uniqueness problems 
are easily solved and calculations can be algorithmically defined. Since we are 
here interested in applications, we shall skip the rigorous geometric framework 
and we refer the reader to [Q Q for a deeper insight into the mathematical 
details of the matter. To stimulate the interest of physically-oriented readers 
we shall keep as much as possible the formal setting to a minimum and we shall 
use a coordinate language. 

Given a vertical vector field X, the variation Sx^ of the Lagrangian L = 
C{j^y) ds can be generally written, through a well known integration by parts 
procedure, as follows: 

5^/:(/y) = e,{f^y)X^ + d>^F\f'^-\j,j'^-'X) (1) 

or, in terms of differential forms, as 

SxL^e{L,X) + dF{L,X) (2) 

Here the m-form 

e{L,X) = e,if''y)X'ds (3) 

and the (m — l)-form 

F{L,X)^F\f''-^y,j''-^X)dsx, ds^=ixds (4) 
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are called, respectively, the Euler-Lagrange form and the Poincare-Cartan form. 
We stress that the Poincare-Cartan form F{L,X) depends linearly on X to- 
gether with its derivatives up to the order k ~ 1 included. [We also remark 
that, for higher order theories the Poincare-Cartan form is not unique since it 
usually depends on the choice of a linear connection on spacetime; see |^. We 
shall not care about this in the following since we shall only deal with General 
Relativity where a canonical Poincare-Cartan form can be chosen.] 

The critical sections a : x^ t-^ y^ = o''^{x) are those which satisfy the Euler- 
Lagrange equations: 

e^if'y)\y^^^^^=0 (5) 

We shall now consider natural theories, namely, theories describing the be- 
havior of geometric objects (such, e.g., tensor fields, tensor densities or linear 
connections) by means of Lagrangians which are covariant with respect to the 
action of diffeomorphisms of spacetime M. In mathematical language natural 
theories are the ones fulfilling the fundam,ental identity 



dy' ^^ '" dyl^...^ 



di'^L) = \ i^^^^y' + ■■■ + 1^. — ^^yi-i..... \ ds (6) 



for each vector field i^ on spacetime; see |1^. Inserting the first variation formula 

(2) (with X replaced by £^y) into the right hand side of the fundamental identity 

(3) we obtain the conservation laws 

d£{L,0 = W{L,0 (7) 

where the {m — l)-form £{L,S,) and the ?n-form W{L,S,) are defined, respec- 
tively, as follows: 

£{L,^) = F{L,£^y)-i^L (8) 

W{L,^) = -e{L,£^y) (9) 

and are called, respectively, the Mother current of L relative to i^ and the work 
current. Notice that, by definition, one has 'W{L,S,)\y^^i^\ = whenever a 
is a solution of field equations, so that >V(L,^) evaluates by (0) the "work" 
performed off-shell, i.e. by sections which are not critical. Accordingly, from 
equation (0) we infer that the differential form £{L,^) is conserved (i.e. closed) 
on-shell. Moreover, in all natural theories the map ^ i-^ £{L,£^) is a linear 
partial differential operator in the coefficients £,^. Hence, the Nother current 
£{L,^), through a (covariant) integration by parts, can be rewritten as 

£(L,e)=£(L,0 + cW(L,0 (10) 

where the (m — l)-form £{L, ^) is called the reduced current and it is vanishing 
on-shell because it is proportional to a combination of field equations; see [^8|, 
|l4| . The (m — 2)-form W(I/, ^) is instead called the superpotential. It depends 
on the fields y and their derivatives up to order 2k — 2 and it is linear in the 



components ^^ and their derivatives up to order fc — 2. Hence, it follows from 
(lOh that the Nother current £ is not only closed but it is also exact on-shell.n 
Since £{L,X) is an {m — l)-form in spacetime, it can be integrated over an 
(m — l)-dimensional region S, namely, a submanifold S of M with a boundary 
dT, C E C M which, in turn, is a compact (m — 2) -dimensional submanifold. 
With this formalism at hands, one could be tempted to define the Nother charge 
Q-s{L, ^, a) along a section a and relative to ^ as follows: 



Q^{L,^,a) = [ £{L,0\ 



V=(j(x) 



£{L,0 ^ +/ Z^(i,6l,=.(.) (11) 

(notice that the Nother charge QeI-^jCi'^) becomes a pure boundary integral 
whenever cr is a solution of field equations). Nevertheless it is well known that 
the definition given in this way does not fit the physical expectation values. 

Example 2.1 The Hilbert Lagrangian 

LH = ^^g'"' R^.uds (12) 

(with K = Stt in geometric units with = 0=1) is a second order Lagrangian 
(fc = 2) on the configuration bundle Lor{M) — > M of Lorcntzian metric on the 
4 dimensional spacetime M . A vertical vector field on the configuration bundle 
reads locally as, X = {5g^^)d/dg^i,u- The variation of the Lagrangian can be 
written as: 



5LH = '^^G^''8g^,ds + d 



2k 



:^g^-5ui,ds. 



(13) 



where 



""ly — r"^"^ ~ ^'{p^Z)a G'^jy — R^u — -j^g^iv R (14) 

Hence, the Eule-Lagrange form and the Poincare-Cartan form are given, re- 
spectively, by 

e(Lff,X) = -i-V^G^''%,ds (15) 



V5, 
2k- 



F{LH.X) = ^^g^^''6u^^,ds^ (16) 



The Nother current, defined as 



£{LH,i) = ^ {.g^^i^cKJ -^Lh) ds^ 



^We stress that the Nother current as well as the superpotential are algorithmically and 
uniquely defined at the bundle level in terms of jet bundle morphisms and they are canonically 
associated to the Lagrangian. Then they are computed along a configuration a and they give 
rise to differential forms on spacetime. 



splits as £{Lhi = £{Lh, + dU{LH, where 

i{LH,0^—G':Cdso. (17) 

K 

U{LH,£.) = ^9'""^peds^, (18) 

Notice that the superpotentialZ-/(L//, S,) is noting but the Koniar superpotential; 
see [p8J. When integrated on spacehke hypersurfaces, it is well-known that the 
Koniar superpotential does not reproduce the expected values for all conserved 
quantities. Indeed it is affected by the so-called "anomalous factor problem" 
when computed along asymptotically flat solutions (see ||2^) or, even worse, 
it gives rise to divergence problems on solutions of Einstein's equations not 
featuring an asymptotically flat behavior (see |3y, |3l|] ) . 

This example strongly suggests that a modification in the definition of Nother 
charges has to be provided. ■ 

Let us consider again a general natural theory. According to definition (^, 
for each vertical vector field X = (^J/') d/dy^ we have: 

5x£{.L,i) = 5xF{L,£^y)-i^{5xL)^^5xF{L,£e,v)~He{L,X)-^dF{L.X) 
= SxF{L, £^y) - £^F{L, X) - i^e{L, X) + di^F{L, X) 
= u{L,X,£e_y)~i^e{L,X) + d{^F{L,X)) (19) 

where the symplectic current lu{L, X,Y) relative to L and calculated for two 
vertical vector fields X and Y is defined as (see p^): 

coiL,X,Y)=dxF{L,Y)-SYF{L,X) (20) 

Equation (|l^) suggests to redefine the variation SxQ of the corrected conserved 
quantity Q as follows 

SxQsiL, ^ [ [Sx£{L, " d{i^F{L, X))] 

'=' / SxSiL, 0+ f [5xU{L, - i^F{L, X)] (21) 
Je Jot. 

Hence, from equation ( [l9| ) we obtain: 

'5xQs(i,e)= j W{L,X,£^y)^i^e{L,X)} (22) 

This prescription is at the basis of the so-called ADM covariant formalism; see 
[O, ^, as well as it is at the basis of a Hamiltonian analysis through a sym- 
plectic viewpoint (see [O, O, ^). Indeed, whenever the (nowhere vanishing) 
vector field ^ is transverse to a Cauchy hypersurface S we define the (variation 
of the) Hamiltonian 7J(L,^,E) simply as follows: 

SxH{L,t^) = SxQj:iL,0 (23) 
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i. e. we identify the variation of the Hamiltonian with the variation of the 
corrected Nother charge Qs relative to £,. First of all notice that definition 
( p3[ ) has the advantage of being independent on divergence terms added to the 
Lagrangian. Indeed if we consider a new Lagrangian L' differing from L by a 
total divergence, i.e. L' = L + d(3, where /3 is a (m — l)-form, we have 

U{L',i)=U{L,S,)+i^P 
^F{L',X)^i^F{L,X) + 1^5x13 

and hence 5xQt.{L' ,S,^) — 6xQsiL,S,)', see definition (|2l|). This means that 
divergence terms added to the Lagrangian just affect the superpotential so that 
this fact would reverberate on the definition of the Hamiltonian if it were naively 
defined through the non corrected Nother charge (111]). Namely, different diver- 
gence terms would imply different definitions of Hamiltonian and, in turn, dif- 
ferent "boundary Hamilton equations". Instead, the boundary term i^F{L^X) 



added in the modified definition ( 21 ) exactly counterbalances the modifications 
in the superpotential and we end up with a unique definition of Hamiltonian 
(unique inside the set of Lagrangians differing from each other only by diver- 
gence terms). This property is clearly welcome. The duty of the Hamiltonian 
should in fact be to dictate, via Hamilton equations, the evolution laws of the 
dynamical fields in such a way that Hamilton equations have to correspond to 
(the space-1-time decomposition of) Euler-Lagrange equations. As well as these 
latter equations are insensitive to Lagrangian divergence terms one would expect 
the same property to hold true also for Hamilton equations. 



Notice also that, according to definitions (21) and (|23D, the (variation of 
the) Hamiltonian consists of a bulk term plus a surface term. It is exactly the 
surface term which plays a key role in testing the viability in the definition of 
the Hamiltonian. Indeed the variation of the Hamiltonian is well-defined if all 
surface terms in 5H vanish. That is, the surface terms in the Hamiltonian, when 
varied, should exactly cancel out surface terms arising in the variation of the 
bulk term. In this way, the variation of the Hamiltonian would be well-defined 
on the full phase space of the dynamical fields rather than just on the restricted 
space of fields satisfying suitable boundary conditions (namely, the boundary 
conditions necessary to cancel out possible surface terms in 5H). Hence, in 
order to give a viable definition of Hamiltonian, we propose here an additional 
boundary modification to the general prescription (pi|). We stress that, when 
dealing with General Relativity, the following modification will turn out to be 
relevant only when the foliation of (a region of) spacetime does not feature 
orthogonal intersections with the timelike boundary. 

Let us thence consider a (portion of a) Cauchy hypersurface S in spacetime. 
Since the Poincare-Cartan form F{L,X) is linear in X and its derivatives, 
owing to definition ( po|) the integral f^ i^iL, X, £^y) of the symplectic form on 
the hypersurface S, by means of an integration by parts procedure, can be 
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rewritten as follows: 

/ uj(L, X, £^v) = f Cj{L, X, £^y) + f t{L, X, £^y) (24) 

i.e. it splits into a bulk term plus a pure boundary term which, for fc > 1, is in 
general non vanishing (and we shall see below that this is indeed the case for 
General Relativity). Accordingly, our proposal is to slightly modify definition 
( pi| ) by pushing the boundary term Jg^ t directly into the definition of the 
Hamiltonian and defining then a new surface-adapted Hamiltonian H{L, ^, E). 
That is, we define the variation 6xH of the new Hamiltonian as follows 

SxHiL,^^:) - J {dx£{L,0-d[^F{L,X) + TmX,£^y)]} (25) 

= / SxSiL, + / [SxU{L, - i^F{L, X) - t{L, X, £^y)] 

This prescription is very close in spirit to the original idea of Regge-Teitelboim 
(see |20|) even though we are now dealing with any natural field theory and 
the prescription is based on a symplectic analysis of Nothcr theorem instead of 
being based on a canonical analysis of the ADM Hamiltonian. It follows in fact 
from (H) and (||) that: 

5xH[L, ^, S) = / Co[L, X, £^y) - / i^e{L, X) (26) 



and equation (26) ensures that 5xH does not contain surface integrals which, 
as already explained, would correspond to Hamiltonian boundary conditions on 
the dynamical fields. 

Notice also that the definition (ESh does not alter the property of invariance 
of 5H with respect to the addition of divergence terms into the Lagrangian. 
Nevertheless the definition given depends in general on the particular hypersur- 
face S we are dealing with. Additional boundary hypotheses have then to be 
imposed in order to avoid the problem. We refer the reader to [p3 where the 
dependence of the symplectic current on the Cauchy hypersurface was carefully 
investigated. We shall here skip the details on the matter and we shall instead 
focus on the definition of energy which is our main concern in this paper. 

We define the variation of the energy E{L, ^) relative to ^ as the on shell value 
of the variation of the Hamiltonian. Hence, the variation SxE(L,^) is obtained 
evaluating ( p5|) along solutions of field equations and considering variations with 
respect to a vertical vector field X which, in turn, is a solution of linearized field 
equations, which in our formalism are equivalent to the equations 

dx£{L,O^0 (27) 



Since under these hypoteses we have i^e{L, X) — 0, from definition ( pSD we 
obtain that SxE{L,^) turns out to be a pure boundary term, i.e. 



SxE{L,0= f 5xU{L,0-^F{L,X)-T{L,X,£^y) 



(28) 
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while equation (|2^) reduces to 

5xE(L,0= f Q{L,X,£^y) (29) 



Notice however that the expressions ( |25| ) and (Eq) define only the variation 
of the Hamiltonian and of the energy, respectively. The problem arises whether 
those expressions can be formally integrated or not. We shall see that, at 
least for General Relativity, this problem admits a positive answer. Roughly 
speaking, the problem of integrability is related to the way in which we keep 
the field variations or, in other words, is related to the way we impose boundary 
conditions (e.g. Dirichlet or Neumann conditions). Notice that in the present 
framework for boundary conditions we actually mean conditions which have to 
be satisfied by the vertical vector fields X on the boundary SS (for example 
X\g-^ = 0). Different classes of vector fields correspond to different sets of 
boundary conditions and they will give rise to different definitions of energy, 
as it is physically expected. Notice that boundary conditions are not imposed, 
a priori, by adding divergences terms to the Lagrangian we started from. As 
already outlined in the Introduction, this would be a rather unpleasant strategy 
mainly because we would have to modify the theory from the very beginning, i.e. 
at the Lagrangian level, accordingly on what kind of energy we are interested in 
(e.g. internal energy, free energy and so on). Indeed we stress again that each 
kind of energy is deeply related to its own boundary conditions (see M, nH, M ) 
and in the Hamiltonian formulations based on the canonical reduction of the 
action functional, boundary conditions are dictated by divergence terms added 
to the Lagrangian and vice versa. Hence, in order to satisfy different boundary 
conditions we would have to add different divergence terms each time. In other 
words, each symplectic control-response boundary structure dictates its own 
action functional. 

On the contrary our Nother-based prescription is insensitive on divergence 
terms so that definition (|2l) is unchenged if we remain inside any class [L] of 
Lagrangians, the elements of which consist of Lagrangians differing from each 
other only by divergence terms. Different physical meanings can be nevertheless 
attributed to the formula (p8|), which is the same for each element of the class [L], 
by considering variations with respect to different classes of vertical vector fields 
X. Namely, boundary conditions play an effective role only at the Hamiltonian 
level. The vector field X (together with its jet prolongation) keeps track of what 
has to be held fixed on the boundary, i.e. it determines the control variables. In 
this way the Lagrangian formulation gets rid of a direct physical meaning which 
is instead restored at the Hamiltonian level. Indeed, it is only in the Hamiltonian 
framework where we read off the symplectic structure of the theory, namely the 
control-response variables and the corresponding energetic content; see [^, |T5[. 
We can completely ignore divergence terms and the corresponding boundary 
conditions in the Lagrangian formulation provided we assign to the Lagrangian 
the only duty to frame the theory in a well-defined geometric framework and 
to provide the equations of motionR. Instead, at the Hamiltonian level, we are 

•^ Notice that, in the geometric approach to Lagrangian theories, the equations of motion are 
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free to fix the boundary conditions in any way we want through, the choice of 
appropriate vector fields. 

The application to General Relativity will help us to clarify in detail how 
the idea works on. 



3 General Relativity 

As an application of the formalism let us now consider again standard General 
Relativity in vacuum with dim M — A (notice however that the formalism goes 
through for each spacetime dimension m > 2. The case m = 2 will be instead 
treated separately in a forthcoming paper Q). 

Let us consider a 4-dimensional region D C AI which is diffeomorphic to 
the product E x K where E is a 3-dimensional closed manifold with boundary 
B — 9S. We denote this diffeomorphism by 

V- : E X 1? — >D (30) 

For simplicity we assume in the present Section that the boundary B — 9E is 
topologically a sphere (even though no relevant change arises when we endow 
E with an internal boundary, i.e. when the manifold E is assumed to be the 
product S^ X I where / is a closed real interval). For any t G M a hypersurface 
Et C I? is induced by i/; according to the rule T,t = ipCS x t) ~ i^tC^) and 
we require it to be a spacelike Cauchy hypersurface. Each map tpt : S — > Ej 
establishes an embedding of E into M. The set of all Et, for t in M, defines a 
foliation of D labelled by the time parameter t. Moreover, each Et intersects 
the boundary dD in a compact 2 dimensional surface Bt which is diffeomorphic 
to B, for all t in M; see Figure 1. The diffeomorphism is established by the map 
tpt ■ B — > Bt- Hence the boundary B = dD is a timelike hypersurface globally 
diffeomorphic to the product manifold B x M. The time evolution field ^ in 
D is defined through the (local) rule ^'^V^i = 1 and, on the boundary B, it is 
tangent to the boundary itself. We shall use Greek letters, ranging from to 3, 
to label indices of tensorial objects on D. Tensors on E are labelled by Latin 
letters a,b, . . . running from 1 to 3. Tensors on B x M are labelled by middle 
Latin letters i,j,--- — 0, 2, 3; while Capital latin letters A,B, . . . = 2,3 are used 
for tensors on B. Notice that tensorial objects T^a'" on M, with contravariant 
indices fijiy . . . tangent to each Ej, can be pulled back on E via the embedding 
maps ipt and they originate tensor fields T^j" on E. Accordingly these tensor 
fields have a double interpretation. When labelled by Greek indices they are 
thought of as tensor fields in spacetime M, whilst when they are labelled by 
Latin indices they are instead thought of as tensor on E. The same double 
notation applies, with the corresponding modifications, for tensor fields which 

mathematically recovered from the variation 5x ^d (cr) of the action functional by postulating 
that critical sections a are the ones extremising the action for any region D and for any vector 
field X with compact support in D. Hence all boundary integrals in SxAo{a) are vanishing 
due to the hypothesis made on X. In this way action functionals differing from each other by 
boundary terms give rise to the same equations of motion. 
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are "tangent" to B and tensor fields onBxlR; and also for tensor fields "tangent" 
to Bt and tensor fields on B (where "tangent" is referred to each controvariant 
index) . 

We denote by u^ the future directed unit normal to Ej and we denote by 
n^ the outward pointing unit normal of Bt in E^. Of course u'^n^j^ =0. 
Notice that we are not requiring that the hypersufaces Et intersect the timelike 
boundary B orthogonally so that, in general, one has £,^nf^\g ^ 0. We also 
denote by n^ the outward pointing unit normal of B in M and with u^ the 
future directed unit normal of Bt \tiB. We have C^ "-/^J^ = u'^ n^|^ =0. The 
evolution vector field can be decomposed as <^'^ = TV u^ + N^ where the shift 
vector iV^ is tangent to the hypersurfaces Ej, i.e. iV^ u^ — (but in general 
-^^"-mIs 7^ 0)- The metrics induced on Ej, S and Bt by the metric g^^ are 
given, respectively, by: 

hyiv = a^jiv + u^^ u„ (31) 

If^i' = gfii' - TT-iJ. riu (32) 

(Notice that hat, Jij and aAB will instead denote the pull back metrics on E, 
B X M and B, respectively.) We shall denote by h, 7 and a, respectively, the 
absolute values of the metric determinants. The metrics (|3l|), (^) and ( ^ 
with an index raised through the contravariant metric g^'^ define the projection 
operators in the corresponding surfaces. We also denote by 

Kf,, = -/i^VaU. (34) 

©M- = -7^V„n, (35) 

ICfj^u = -a'^DaU^ (36) 

respectively, the extrinsic curvatures of Ej in M, of K in M and of Bt in E^. 
The symbols D and T) denote, respectively, the (metric) covariant derivative on 
Et compatible with /i^,y and the covariant derivative on B compatible with the 
metric 7^^. We shaU denote by if = K^^h'"', 6 = Q^^^'"' and /C = K.f.^a'"' 
the traces of the appropriate extrinsic curvatures. The momentum P'^'^ of the 
hypersurface Ej is defined as: 

P^"' = — {Kh>"' - Rt"") (37) 

while for the momentum H^'' of the hypersurface B we set: 

n'"' = - ^ (67^'' - e^") (38) 

From a physical viewpoint the surface Bt can be thought of as an istantaneous 
configuration of a set of observers at the time t and their history is represented 
by the timelike boundary B; see [H. The unboosted, or unbarred, observers are 
the ones evolving with four velocity w^ and hence orthogonally to the foliation 
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Et. They are at rest with respect to the S fohation. They view m'^ and n^ as the 
unit normals to Bt . We instead refer to the boosted, or barred, observers to the 
ones which are on the same surface Bt at the istant of time t but evolving with 
velocity ^'^. They are comoving with B. Since ^'^ € TB then boosted observers 
view u^ and n^ as the unit normal to Bt. Barred and unbarred observers are 
related by the boost relations: 



u 



jUf^ + -fv n^ (39) 



P f "^^l I I ^ '"p, 



n 



M 



7 n^ + 7W Uf, (40) 



where 

7W = u^n'^ = —u^n^ = sinhfl (41) 

measures the non-orthogonality of the intersections of the leaves S( with the 
timelike boundary B. The parameter 9 is usually called the velocity parameter. 
Moreover we set 7 = (1 — u^)"^/^ where 

^^JJ^^N^ (42) 

is the boost velocity between the observers. It measures the velocity in the 
direction of n^ of an object with four velocity ^^ as measured by the observer 
with velocity u^, namely the velocity in the n^ direction of boosted observers 
with respect to unboosted observers. 

The vector field 1^'^, as already said, can be decomposed as ^^ = N u^ -\- N^. 
On the boundary B it can be also decomposed as ^'^ = N u^ + N^ where 
N = N/^ and 7V^ — a'^N^ . From now on we shall use the notation that a 
bar over a quantity will denote that it is defined with respect to the normals 
u^ and n^ , so that, e.g., IC^i, would correspond to the extrinsic curvature of Bt 
embedded in a spacelike hypersurface with unit normal u^ in Bt (or, roughly 
speaking, it is the extrinsic curvature of Bt with respect to the normal n^). 
Moreover, the following identities hold true: 

^=NVh, Vi = NV^ (43) 

3.1 Hamiltonian 

Let us now consider the Hilbert Lagrangian Lh- According to the general 



prescription (25) the variation of the Hamiltonian is defined as follows: 



dxH{LH,C,^t) = / Sx£iLH,0 (44) 



/ [6xUk{Lh, - hF{Lh, X) - t{Lh,X, £^g)] 

JBt 



where (see example 2.1) we set: 



2k- 



F{LH,X) = ^^g^^''5u^^,ds^ (45) 
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£(Lh,S.) ^^^GZ rrfsa = {NU + iV"Ha} <r X (46) 

2k- 



^i^(^ff'0 = &'^''Vpe^ds^. (47) 



and Ji and Tia in ( |46| ) are called, respectively, the Hamiltonian constraint and 
the diffeomorphism constraints. 

Now, our next task will be to calculate the term t{Lh, X, £^ g) and to 
rewrite ( [44| ) in terms of quantities adapted to the spacetime foliation. In this 
way it will be possible to establish whether formula ( [44| ) is a viable definition 
for the Hamiltonian. Moreover a correspondence with other definitions exist- 
ing in literature will be possible. From now on we shall adopt the notational 
choices made in ||] and we shall make repeated use of useful formulae there 
demonstrated. Let us start with the relations 

ff"''K/3 "m = 2<5i^ + K'^'Sh^, + D^ih^Su'') (48) 

9"^^<p ^M = 2Je + e^^'^Sj^, + V^irjnn (49) 

(the proof of which can be found in the aforementioned paper g]). From ([45|), 
( [isl ) and Stokes' theorem it then follows that 

2k Jst 2k Jg^ 



[ h^JP^'^'d^x -^ I n^,8u^^fSd^x 
Jst 2k Jst 



(50) 



(notice that Jj, ^ /^ ds^ = — J^ f^ Up V^ d-^a; for any 3-form f — ft^ ds^ 
and /j, Dfj^X^^ \/hd^x — J^ X^^ n^^/a d'^x for a vector field X'^). Let us now 
consider the integral of the symplectic current ( po| ) on the generic leaf Et of the 
foliation. From ( pO[) we obtain: 



/ uj{LH,X,£^g) = u}{LH,X,£^g)+ T{LH,X,£^g) 



(51) 



where 



c^(Lff , X, £^g) = {(i;^ V) '5^'"' - (i^e^^'^) 5h^,} d^ X (52) 

t{Lh,X, £^g) = ^ {i:^(\/5^n^5w'') - S{y/^n^£^u^')} d^ x (53) 

Notice that, having not imposed the hypothesis of orthogonal intersections, 
the 2-form t(Lh, X, £^g) does not vanish at the boundary Bt since n^ in 
general is not surface-forming (instead, in presence of orthogonal boundaries, 
Ufj, is surface-forming and hence dn^, is proportional to n^ so that n^ Su'^ = 
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—u^Sn^i, = 0). Hence formula (^ actually constitutes a generalization of pre- 
vious definitions found in literature and based on Nother theorem; see [12, O, 

The direct calculation of the expression (^ turns out to be rather cum- 
bersome. Fortunately many results can be found elsewhere and we shall only 
reproduce the final results. First of all the bulk term in (El) has been already 
computed in 0, 0] even though it was originally calculated with respect to un- 
barred quantities. In terms of barred quantities it reads as follows: 

/ 5x£{LH.i,^t) = / {5Nn + 5N'^no. + [ho.i3\5P'''^ -[P'''']5h^p]d'x 

Bt I ^ 

+ -^[£^{V^)5{e)-5{^){£^6)\\ (54) 

where 

e = -£ (55) 

~jo. = -^a^^M^"'u, (56) 

V7 

r'^ = -[(n''a^)a"'^-£a"''+,^"'^] {a^=u''V,u^) (57) 

are, respectively, the surface quasilocal energy, the surface angular momentum 
and the surface stress tensor, as viewed by the boosted observers. They are all 
tensors on the boundary surface Bt and, after integration on the surface, they 
can be used to describe the stress energy-momentum content of the gravitational 
field inside Bt; see [|l|, ||, ^ ^ and reference quoted therein. Instead [hap] and 
[P"'''^] in (H) are a shortcut for 

2k N 
[hcfs] = -j^{2P^p~KpP) + D^N0 + DpNc. (58) 

2k 2 v/i 

_27Vk . _ l/2PP'^'') + 

+ yjl(^D"D'^N - h"'^h'"'Df,D^N) + 
2k 

+Df,{P'^'^N^')-Di,N°'P^^'-Df,N'^P°'f' (59) 

having denoted by TZaf3 the Ricci tensor of the 3-metric hap. 

Let us now consider the surface terms into (H) to which we shall refer as 
SxH JLH,£,,'St)- We shall carry out the calculation referring to the set of 
barred observers (the ones comoving with B) instead of referring to unbarred 



18 



observers (at rest with respect to the S fohation). We shall also make repeated 
use of the property: 

U = f ^d?xU^°''^^n^up = I ^(fxU^°''^'^ncUp (60) 

Bt ^ Bt '^ Bt 

which holds true for any 2-form U = 1/2 JjjU°'^ dsap- 

The first term in the surface integral (44), after a bit of algebra (see 11^, 0|) 
can be written as follows 



Uk{Lh,0 



1 
2k 



d?x5[2^/^u''Q%a\ 



(61) 



(62) 



+ -- / (fx5 \^/aua£en°'^ 
For the second term in 5xH ,ALh,S., St) we obtain instead: 

- f i^F{Lh,S,) = ^ ( ^ffxN[25Q + Q^'''5^^,] 

J Bt ^'* -I Bt 

The third surface term /gT has already been written down in (|5^), i.e. 

-( T{LH,X,£^g)^^ f {~£^{y^n^6u'') + S{y^n^£^u'')}d'^x (63) 

JBt ^'* JBt 

The value of SxH^^^^{LH,^,^t) results by the sum of (^, (^) and ( |63|) . In 
order to gain some simplification we shall make use of the relations (sec H ) : 

6^11^ = ICf_tu + {n°'aa)Uf,u^ + 2(Tf-^u^-jKa0n'^ (64) 

S%, = -i2/N)u^uJN - (2/7V)a„(^u,),5iV" + af^aj^a^^ (65) 
n^du'' + u^Sn^" = -269 (66) 



where 

After a little of algebra we reach the final result: 



(67) 



6xHALH,^,^t) 



d?x I N5{^1) - N^'SiV^j^) + ^^s'^^Sa^p 



+ - [£dV^) S{e) - S{V^) {£^6)] (68) 



Finally, the sum of the Hamiltonian bulk term (BJ) and the Hamiltonian surface 
term (|6^) gives the result 



Sx 



H{Lh, e, St) = / {SNH + SN'^Ho. + [Kt}\5P''^ - [P^^^l^J/ia/s} d?x (69) 
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We stress that the boundary terms arising in the variation of the bulk term in 
(p4[) exactly cancel out the variation of the boundary term (pSh . 



In order to explicitly write down the Hamilton equations (Eq) wc still have 
to calculate the term 

/ t^e(LH,X) = -^ f ^G^^'^Sg^A'dsp (70) 

Taking into account the relations 






we have 



/ i^e{LH,X)= (71) 

Taking into account (p2|), (p9h and (|7l|), the Hamilton equations (Eq), i.e.: 






(72) 



read as follows: 

Vh 



K 



G^^u^u^ =n = (73) 



-— G^^V^-^Ha^O (74) 

K 

£^h^_,u = {h^,v\ (75) 

^G'^^h'^h^ = £^P^"' - [P^''] = (76) 

Notice that equations (|73|), (|74| ) and (JT^) correspond to Einstein's equation in 
the (3 + 1) formalism while ( |75| ) is nothing but the definition of the momentum 
P^^'^, i.e. the Legendre transformation. 

Hence the prescription (B3) originates the correct evolution equations in 
the phase space {hab, P°-^) of General Relativity and it does not give rise to 
additional boundary equations; see g, |3^. Q For this reason the definition ( p| ) 

^Notice, however, that Hainihon boundary equations would indeed arise if the correction 
term J^ r in formula ( 14J) were missing. 
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can be truly considered as a genuine generalization of the Regge-Teitelboim 
Hamiltonian formulation pO| to arbitrary spacetime dimensions and to solutions 
with arbitrary asymptotic behavior (and not necessarily stationary). 

Notice also the close formal analogy of expression (f2|) with the definition 
of Hamiltonian in Classical Mechanics. In Classical Mechanics one defines the 
generalized Hamiltonian through the usual rule H{q,p)\ , ■^ = p{q, q) q—L{q, q) 
so that 

6H = qSp - p6q +{p- dL/dq) 6q (77) 

In trying to establish an analogy between (^) and definition ^73) the term 
/j, uj{Lh, C^t) in ( [72[ ) can be seen as the qSp — p5q term - see (|52|) - while the 
term J^, i^e(L/f , X) clearly corresponds to Lagrange equations (p — dL/dq) Sq 
- see @). 

The question arises whether a Hamiltonian H exists such that its variation 
5H equals the expression (£^). The answer to the problem is given by: 

H{LH,^,^t)-Ho{LH,^,^t)= f {iVW + iV"H„}d'x (78) 



Bt 



/^{Ne-N°'ja}(fx+- I £^{y/^)e(fx 



where Ho{LH,S,,^t) turns out to be a constant of integration in the left hand 
side of (|6^). [It will play a fundamental role when dealing with the definition of 
energy since it will play the role of a reference term which defines the zero level 
for the energy itself]. Indeed, the variation of ( |78| ) turns out to be: 

SxH{LH,^,^t) = / {6NH + SN''no. + Kp]SP"^'-[P"'']Sh^p}d^x 

-/ iu^''S%,-^£^{e6V^)\d^x (79) 



and it gives rise to expression ( p9D provided that, during the variation, the 
boundary metric 7^^ is kept fixed. [We also remark that it is indeed the Dirichlet 
condition (57 |g = together with the requirement that ^ be a boundary Killing 
vector field that ensures the definition of 6H to be independent on the choice 
of a particular Cauchy hypersurface] . 

Moreover, expression ( |79|) can be used to calculate the time rate of change 
of the Hamiltonian. To this end we have to replace the variations S of the 
dynamical variables in (^9|) with their Lie dragging £5 along the evolution vector 
field ^. Assuming that Einstein's equations are satisfied we obtain: 

£iH{LHA,^t)^- f (n''-'£i%,~^£^{e£^V^)\<fx (80) 

which is vanishing if £(^ 7^,^ — 0, i.e. if the vector field ^ is a boundary Killing 
vector field for the metric 7 or, in other words, if the gravitational system is 
isolated from the outside; see |3, 0, 0]. 

21 



3.2 Energy in General Relativity 

Let us now consider the definition of energy. It could be obtained simply by 
evaluating ( |78| ) on-shell. In order to better understand the role played by the 
reference term Hq and its relations with boundary conditions let us consider 
again the on-shell variation of the Hamiltonian. We specialize definition ( p8| ) 
for the Hilbert Lagrangian Lh so that 6xE{g) = SxH {LH,(,,^t) and from 
( pSi ) one has: 

SxE{g) = J d^xfNS{V^e)~N"6{V^jo^) + ^^s"^Sa^f, 

+ i[£^{^)S{e)~5iV^)i£^0)]^ (81) 

Notice that SxE{g) is a pure surface integral on Bf. Nevertheless we stress that 
definition ( pi| ) is not a quasilocal energy relative to the surface Bt- Indeed, 
let us now consider a hypersurface St with an external boundary Bt and an 
internal boundary Ct, i.e. 9Sf — Bt(J Ct- Taking expression (|5^ ) into account, 
equation (129) now reads as: 

SxE{g)\Q^^ = SxE{g)\s^-SxE{g)\c^ 

[(i:^ V) ^P^" - {£iPn Sh^.] d^'x (82) 

St 

For stationary solutions the right hand side vanishes together with the correction 
term t{Lh, X^ ^), so that SxE has the same numerical value when computed on 
different surfaces belonging to the same homology class. In particular, if spatial 
infinity is homologic to the surface Bt, expression (|8l| ) corresponds to the total 
energy or, better, it corresponds to the total mass. Applications to a number 
of specific solutions, including e.g. Kerr-Newman, BTZ, Taub-bolt solutions, 
can be found in jlj, pO^, |3l| . In all the examples so far analysed it has been 
shown that expression (|8l|) (with r = 0) always gives rise to the expected value, 
namely, the total mass of the solution. 

As previously remarked an expression of the quasilocal energy can be instead 
obtained only when boundary conditions are imposed on the surface Bt- We 
define the quasilocal internal energy E(g, Bt) contained in the surface Bt to be 
the value obtained from ( pl| ) by imposing Dirichlet metric boundary conditions. 
Namely, we consider a one parameter family of solutions of field equations which 
admit the same boundary metric 7^1/ and we consider a vector field X tangent 
to this family, i.e. 

S%.\b = (83) 

Since ^7p,y|g = impHes SN\j^ = 0, SN"\^ = and (5(T^,y|g = 0, formula ( ^ 
can be explicitly integrated: 

E{g,Bt)-Eo{go,Bt)^ J d^x fv^ {N e- N^'jo.) + ^£^{V^)e\ (84) 



22 



We stress that we are integrating along a family of solutions satisfying the 
boundary condition ( |83| ) and thence the term EQ{gQ,Bt) corresponds to the 
quasilocal energy relative to a solution go inside the family itself, which then be- 
comes the reference point (or zero level) for the energy. Owing to the properties 
N = No\„, N" = Nq\ and cr^j, = croAii'le: expression (^ can be equivalently 
rewritten as: 

E{g, Bt) =J d^xl/^ [N{e - eo) - iV"(J„ - joa)] + ^£dV^) {0 ~ do)] (85) 

where the subscript clearly refers to the background solution gQ. In other 
words, the background terms can be considered as resulting from the embed- 
ding of {Bt^a^v) into a reference spacetime and by requiring that the surface 
evolution is the same in M and in the reference spacetime, i.e. by requiring 
that the boundary lapses N and Nq coincide, the boundary shifts TV" and Ng 
are equal as well as the time rates of change £^{^Ja) of the metric a in the 
two spacetimes are the same; see H|. Roughly speaking, expression ( pq ) de- 
fines the quasilocal energy of the solution g relative to a solution g^ satisfying 
the same Dirichlet conditions on the boundary. From a symplectic viewpoint 
we are selecting the six components of the boundary metric 7^^ as the control 
parameters of the gravitational system, while the response variables are the 
relative quasilocal surface energy e — eg, the relative quasilocal surface angular 
momentum ja — joa and the relative velocity parameter 9 — 9q. 

We stress that formula ( p5| ) differs from the expression of quasilocal internal 
energy as given in Q only for the addition of the term £^{'\/a) {9 — 9q). [Indeed, 
in y the term 9 £^{y/a) was considered as a pg boundary term and thereby 
it was not included in the definition itself of the Hamiltonian] . Hence formula 
(pq) perfectly agrees with the expression given in those papers when the control 
parameter ^/a is constant in time, i.e. £^{\/a) = (or when the intersection 
parameter 9 equals the intersection parameter 9q in the reference spacetime). 
Otherwise expression ( p5| ) contains an additional contribution to the energy 
coming from the accretion in time of the surface. However, we stress that the 
quasilocal energy contained in Bt depends not only on Bt itself but also on 
its history, geometrically described by the world tube B. Notice that we can 
always make use of the condition £^{-^) = to assign a world tube B to Bt 
or, equivalently, to assign an evolution vector field ^ to _Bt; see [H. In practice, 
however, one always assumes that the evolution vector field is a Killing vector 
field on the boundary. Indeed, this requirement guarantees the conservaton of 
energy in time, where the time is defined by the flow of the vector field. Provided 
that these remarks are taken into account, the term £^{^/a) {9 — 9o) can be set 



equal to zero so that formula (85) agrees with the one calculated in M. 

Another expression for the energy of the gravitational field in the case of a 
non-orthogonal foliation of spacetime was previuosly calculated in |^. The 
calculation there developed was carried out with respect to the foliation of 
spacetime as a whole rather then with respect to the foliation of the boundary, 
as instead explicitly done in this Section. 
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The expression there presented, rewritten according to our notation, reads 
as follows: 



Eh nig) ~ EnHigo) 



d x\fa 



Bt 



NK. - n.N^ 









The background solution go was then chosen in B by requiring that g and qq 
induce the same 3-metric on the boundary B and the same intersection angle 9 
on Bt- Hence: 



EHHig) = - f d^xV^N[IC-ICo] 

^ JBt 



(86) 



-2 / d^Xy/^ 
JBt 

Taking into account the relation 



n^iN^ 



Vh 



nof^cNoi, 






(87) 



together with the splitting (see g) 

Kf,,, = If,,, + (u ■ b)nf, n^ + 2 a^^ n^) Ka^p n^ (b^ = n^V^n'") 

(where If^i, = —a'j^a^VaUp and 6^ = n'^V^n'^), expression ( pq ) can be further 
rewritten as: 

EHHig) ^ J d^x^i^[iIC-JCo) + v{l-lo)]-N^a'il{jc,~jo<,)\ (89) 

(with I = a^^" If^u and Iq = a^" /o^i/ = cr'^" kiiu)- 

In order to establish a correspondence between expression ( pq ) and the defi- 
nition (|89|), we first have to rewrite (^) in terms of "unboosted" quantities. To 
reach this goal we shall make use of the boost relations (see again ||^): 



N ^ Nh 

IC ^ J K- + jvl 
1. 



(90) 
(91) 

(92) 



Jq = ja ^aO 

K 

(where Jq, — ^cr^ Kp^j^n^^ = — jr'^a^j.P^'^ny) together with the relation 

£^{^/^) = ^/^[-Nl-NvK, + 'i:)o,{a'^N^)\ (93) 

where S)q. denotes the (metric) covariant derivative on the surface. After a 
straightforward calculation it is easy to show that expression (B5|) exactly agrees 
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with (p% provided the same matching conditions are imposed between g and 
the background solution go- 

For the sake of completeness we recall that a formula for the energy of the 
gravitational field in the case of a non-orthogonal foliation of spacetime can be 
also found in H . It was there obtained through a new method of variation of the 
gravitational Lagrangian and through a Legendre transformation tecnique. We 
end this section by noticing that expressions (^ll) and (p2|) together reproduce 
the "homogenous" formula (80) of Kijowski's paper |^ which is at the basis of 
his formalism. 

In this way a (conditioned) correspondence among the various definitions of 
energy iQ, |[, 0, |§1 and [g| (based, respectively, on the canonical analysis, on 
Hamilton-Jacobi analysis and on symplectic analysis) and the Nother approach 
is established. 

3.3 Entropy 

In the previous section we have shown that formula (|8|) leads to the definition 
of quasilocal internal energy provided that Dirichlet boundary conditions are 
imposed on the metric field. Formula (|8l| ) can be also used in order to describe 
the first law of black hole thermodynamics. Indeed, the direct evaluation of ( pi|) 
on the horizon of a black hole solution and the boundary conditions dictated by 
the geometry of the black hole horizon lead to the formulation of the first law 
of black hole thermodynamics. 

Indeed, let us specialize expression ( pl|) to an axisymmetric stationary black 
hole solution (7^,^ of Einstein (vacuum) field equations. On each leaf Et let us 
choose spatial coordinates x* co-rotating with the horizon. Hence the vector 
field £, coincides with the null vector field generating the horizon. We recall 
that, with this choice of coordinates, the metric near the horizon satisfies the 
regularity condition: 

n^df^N = kh {kh = surface gravity) (94) 

which ensues from the necessity of removing the conical singularity in the com- 
plex metric obtained from the Lorentzian metric by analytically continuing time 
to imaginary values (see [^, Q, |^). We claim that formula ( ^ ) already en- 
closes all the information about the first law, in the same way it was originally 
formulated in B , but now generalized to non-orthogonal intersections (see also 
[Eol). To realize this, first of all notice that the integrand in the right hand side 



of formula (81) is a closed 2-form as ^ is a Killing vector field of the solution 
g (and hence the right hand side of (|8^) is vanishing). By Stokes' theorem, we 
have that 6xE{g)\g — SxE{g)\^ if Ct is any 2-surface homologous to Bt, i.e. 
if Bt — Ct is the boundary of a 3-dimensional region E^. We can choose, for 
example, Ct = A, where A denotes the (cross section of the) black hole horizon. 
In that case the lapse and the shift are vanishing on the horizon and, taking the 
boundary condition (M) into account the only term surviving in the right hand 
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side of (|81|) is (see |0[): 

5xE{g)\^ = ^5 U d'xV^] = ^SAh (95) 

where A^ denotes the area of the horizon. In order to attribute a thermo- 
dynamical meaning to the dynamics of black holes, in accordance with O] we 
define the following intensive variables: the inverse temperature j3 of the surface, 
the proper angular velocity to and the surface pressure p^^ . They are defined, 
respectively, as follows: 

/3 = ^N 
f3u; = —N'^ (96) 

KH 

' 2 ^ KH 2 



As one can read off from ( |8l| ) they correspond to 2ti / kh times the quantities 
conjugated, respectively, to \/^e, \fo ja and Oaj^, i. e. conjugated to the exten- 
sive variables (however, we stress that the physical interpretation of quasilocal 
quantities in terms of thermodynamical variables comes from a path integral 
formulation of gravity - see |g^, g| - and it is then out of the scopes of our 
discussion). Equating the left hand sides of equations (|8l|) with ( p5| ) and tak- 
ing the property of stationarity into account we then obtain the first law of 
black hole thermodynamics for spatially bounded system as formulated in 
but generalized now to non-orthogonal intersections: 

Sx |/3<5(e V^) - ^^Ki^) + fi^v''^^^^ = ^ (97) 



(in geometric units G = c = 1 where k — Stt). As already remarked in [^1, the 
equation ( p7| ) resembles the thermodynamical law T5S = 5E — ujSJ + pSV of 
thermodynamical systems provided we define f3 = 1/T and S — Ah/ 4: (and it 
reduces exactly to it in static, spherically symmetric systems when /?, uj and p*-' 
can be pulled out of the integral) . 

We point out that the equality in formula (|9^) is a trivial consequence of 
homological properties, namely of Stokes' theorem. For this reason the formula 
can be generalized to solutions other than black hole ones and, of course, it does 
not rely on horizons. Additional boundaries (i.e. additional singularities other 
than a single black hole one) are responsible of additional contributions to the 
entropy. In this way, the homological nature of entropy becomes evident; see 

[ID, @' H- 

We also stress that formula (pn) is also well-suited to describe the thermo- 
dynamics of isolated horizons; see [Ell. Indeed evaluating (pW on a isolated 
horizon and taking into account the boundary conditions dictated by the defini- 
tion itself of isolated horizons, the first principle is recovered. We shall analyse 
this matter in the forthcoming paper B3]. 
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Figure 1: non-orthogonal foliation of a spacetime region D 



